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Abstract: We calculate both analytically and numerically the evolution of a highly rel- 
ativistic fireball through the stages of free expansion and coasting, and determine the 
dependence of the thermodynamic and radiation variables in the co moving and laboratory 
frames. Magnetic fields may have been important in the original impulsive event. We 
discuss their effect on the fireball dynamics, and consider also their effects on the radiation 
emitted when the fireball runs into an external medium and is decelerated. The inverse 
synchro-Compton mechanism can then yield high radiative efficiency in the reverse shock, 
producing a burst of non-thermal radiation mainly in the MeV to GeV range whose total 
energy and duration agree with those of typical cosmic gamma ray bursts. 

1. Introduction 

Gamma-ray burst sources (GRB) have long been suspected to originate from the 
sudden release of energy in small regions of space, where the initial energy density and 
characteristic photon energy is so large that an opaque fireball forms (e.g. Cavallo and 
Rees, 1978, Paczyhski , 1986, Goodman, 1986, Shemi and Piran, 1990). If these events 
were located at distances comparable to (or larger than) typical galactic scales, the 
fireballs would necessarily remain optically thick out to radii of at least ~ 10 9 cm (e.g. 
Zdiarzki, 1982, Imamura and Epstein, 1987), depending on the amount of normal electrons 
and baryons mixed in with the pairs. If the fireball originated in a region smaller than this, 
the radiation pressure on the optically thick fireball would cause it to expand; the evolution 
of this fireball was initially thought to lead, when it became optically thin, to the observed 
gamma-ray bursts. This model has recently received increased attention because the spatial 
distribution of GRBs revealed by the BATSE experiment on the Compton Observatory 
strongly suggested an extended galactic halo or a cosmological origin (Meegan, et ai, 1991, 
Fishman, 1992, Paczyhski , 1992). 

In its simplest version, the fireball model failed to account for the duration and 
time-structure of the observed bursts: it predicts a very short burst, emitted when the 
expanding fireball becomes optically thin (Goodman, 1986), and a quasi-thermal, soft 
7-ray spectrum (Goodman, 1986, Paczyhski , 1986). Another problem, emphasized by 
Paczyhski , 1990, was the possibility that the 7-rays could easily be degraded to even 
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lower energies by adiabatic cooling due to scattering on electrons associated with baryons 
polluting the pair flow. In fact, in the latter case, most of the fireball energy gets converted 
into bulk kinetic energy of the polluting baryons, while the photon burst at thinning has 
much less energy than initially produced. These problems are removed when one considers 
anisotropic scenarios where high Lorentz factor fireballs can escape the inevitable baryon- 
polluted slow wind that must accompany the initial energy deposition (e.g. Meszaros and 
Rees, 1992a, 1992b). In such cases, the main part of the observed gamma-ray radiation 
occurs when the bulk kinetic energy of the baryons (carrying essentially the full energy 
of the initial fireball) is re-randomized and radiated away in the blast wave being pushed 
ahead of the relativistically expanding baryons, and in the reverse shock that propagates 
inwards into the baryonic fireball gas as the latter is decelerated by the external medium 
(Rees and Meszaros, 1992, Meszaros and Rees, 1993). This is a very generic mechanism, 
which operates in almost any scenario for the original energy deposition, as long as it 
makes a relatively clean, high entropy initial fireball. Moreover, the process depends on 
the external environment, thereby allowing the possibility that even a standardised type 
of fireball could create bursts with a variety of complex time-profiles. 

The initial development of the fireball determines the final properties of the blast 
waves, and thus of the bursts. In particular the final bulk Lorentz factor, together with 
the external density, determines the blast wave burst duration, while the dynamics of the 
fireball as a function of the baryon loading determines the relative amount of energy in 
bulk kinetic form (which gets radiated in the blast wave and reverse shock) and in pair and 
radiation form (which escapes when the fireball thins, usually before the blast wave burst). 
Since the relativistic dynamics of the fireball expansion determines the entire energetics as 
well as the temporal characteristics of the burst, a detailed calculation is important. 

In the present paper we discuss the dynamics of the relativistic fireball expansion 
from the acceleration stage through the coasting phase, calculating both numerically and 
analytically the thermodynamic and radiation variables of the flow through the transition 
to optical thinness and the saturation of the bulk velocity. We discuss the kinematics both 
in the co moving and laboratory frames, and indicate the scaling of the behavior with the 
various parameters of the problem, extending the treatment to the anisotropic (jet) case. 
We also consider the role of magnetic fields in the dynamics of the evolution, including 
the case of magnetically dominated fireballs, and investigate the effect of the magnetic 
field on the radiative efficiency of the reverse shock that arises when the coasting fireball 
is decelerated in an external medium. 

2. Numerical Treatment of a Highly Relativistic Expanding Gas 

2.1. Computational Approach 

Because of the presence of a natural length in the problem (the initial value of 
the radius r ) the gas dynamic equations do not allow an exact similarity solution valid 
at all radii. For this reason, a numerical solution is the only way to follow exactly the 
development of the gas, although approximate analytic solutions are possible in the initial 
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acceleration stage and in the later coasting phase (see §3). Previous calculations of the free 
expansion of a relativistic gas in spherical symmetry (e.g. Vitello and Salvati, 1976, using 
a characteristics method) have followed the evolution over one and a half decades in the 
expansion factor. However, expansion over such a limited range is usually not sufficient 
for the bulk velocity to reach its ultimate saturation value, for the cases of large values of 
the initial radiation to rest mass energy ratio rj considered in many problems, such as the 
GRB fireball problem, 

rj = ~ 10 3 £ 5 i(M o /0.6 x lO^M©)" 1 . (2.1) 

Here E ~ 10 51 i?5i ergs is the order of magnitude of the initial photon energy expected 
from the liberation of about a solar mass equivalent of gravitational energy (GM & 2 / ~ 
eM Q c 2 ~ 10 54 ergs, where e ^ 1, R n ~ 10 6 cm, most of which goes into neutrinos and 
gravitational waves), while M D is the total mass of polluting baryons that get mixed in the 
photon-pair fireball of energy E a . 

The numerical treatment of the free expansion of a highly relativistic gas with spher- 
ical symmetry calls for the use of Lagrangian hydrodynamics because of the vast range of 
physical scales present in the problem, of the order of ten or more decades. We have devel- 
oped a Lagrangian code for relativistic fluid flows that uses a second order Runge-Kutta 
integrator (Benz 1984) with adaptive time step. This method has the advantage that a 
Courant (Courant, Friedrichs & Lewy 1928) condition does not limit the size of the time 
steps, an essential ingredient to achieve short and large time scales. Relativistic artificial 
viscous stress was implemented in the code and calibrated using relativistic shock tube 
problems. This was used to verify the accuracy of the same code when artificial viscosity 
was switched off, which turns to be quite good (and significantly faster) in the cases of free 
expansion treated here. The reason why this is an excellent approximation is that shocks 
do not normally occur in such free expansion problems. All the runs described below are 
therefore for zero artificial viscosity. Typically, simulations consisted of 200 Lagrangian 
grid points (mass shells). Runs with different number of mass shells were done to check 
the convergence of our results. We monitored the ability of the code to reproduce the 
analytic solutions of the plane symmetric rarefaction wave traveling inwards (Vitello & 
Salvati 1976). We found that for 200 grid points the maximum absolute error (< 3%) as 
well as the maximum cumulative error (< 2%) appeared in the velocity. 

2.2 Numerical Results 

The initial gas configuration adopted here, for simplicity, is that the gas at t = is 
uniformly distributed within r < r Q , having a total rest mass M , and that at t = an 
amount of radiation energy E ^> M c 2 is deposited uniformly throughout the spherical 
gas. The gas is optically thick to its own electrons, and the resulting photon-pair-electron- 
baryon fluid is essentially an isentropic fluid which can be modeled as a gas of adiabatic 
index 7 a = 4/3, as long as the photons are coupled to the matter by radiation drag. The 
same applies to the case where the energy is mainly in the form of magnetic fields (see §4). 
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After decoupling, the adiabatic index becomes 7 a = 5/3. We have carried out calculations 
for values of rj ranging from close to unity up to 10 10 . As a specific example, we discuss 
the case of r\ = 10 4 below, which is characteristic of the relativistic behavior. 

In the course of the expansion most of the mass, as seen by a laboratory observer, is 
concentrated in a thin shell near the leading edge of the expanding gas (c.f. also Vitello 
and Salvati, 1976). This is seen in Fig. 1, where the distribution of the lab density is 
plotted as a function of the (normalized) lab radius for rj = 10 4 . The curves are plotted at 
the times (from top to bottom) when the lab frame expansion factor r(p mx )/r (which is 
proportional to the lab frame time) is equal to 0.2, 0.4, 0.6, 0.8, 1.0 times rj. 

As the gas expands, the mass shells acquire a bulk velocity which initially increases 
in time and eventually saturates. The distribution of this bulk Lorentz factor over the 
various mass shells is shown in Fig. 2, for a fireball of r\ = 10 4 at the lab times (increasing 
from bottom to top) when r(p max ) /r a = rj 1 / 2 , rj, rj 3 / 2 , rj 2 , rj 5 / 2 , rf. One sees that the leading 
edge reaches a value of Y ~ rj, and most of the mass reaches one third of that value, at 
the laboratory time for which the expansion factor r(p mx )/r ~ rj. After that the average 
bulk Lorentz factor of most of the matter saturates to the value V ~ rj, while the leading 
edge remains within a factor 2 of that value, and the trailing edge of the inner 10% of the 
matter tapers off fairly steeply to zero. 

The average bulk Lorentz factor (or what is nearly the same, T(p rnax ) where the 
density is maximal) grows initially linearly with radius, and saturates to a value r ~ rj 
after a (lab) expansion factor equal to rj. This is seen in Fig. 3, where T(p max ) is plotted 
against the expansion factor r(p max )/r . 

Within the thin shell containing most of the mass, called the mass envelope shell, the 
bulk of the mass (say the inner 80% away from the outer and inner edges) is distributed 
quite uniformly, even after the average bulk Lorentz factor has saturated (for r(p max )/r ^ 
rj). Only later, particularly after r(p max )/r ^ rj 2 , a slight asymmetry favoring the leading 
edge of the envelope shell starts to become apparent, but even then the matter within 
the shell can be considered uniform to a good approximation. This is shown in Fig. 4a, 
where the lab frame density Pl is plotted against the Lagrangian mass coordinate M/M t 
(where M t = M a is the total rest mass, which is mostly within the shell) . In figure 4b we 
show the same plot for the co moving mass density p c . The latter is even more uniform 
than the laboratory frame density. Thus, a co moving observer would find that the world 
around it is essentially isotropic and homogeneous, at least far from the edges, but with a 
density that drops in time, as seen from the decreasing value of the average density. The 
various curves, from top to bottom, are for (increasing) times when the expansion factor 
r(pmx)/r is equal to rj 1 / 2 , rj, rj 3 / 2 , rj 2 , rj 5 / 2 , rj 3 . 

The width of the mass envelope shell in the lab frame is initially nearly constant an 
equal to the initial radius, Ar/r Q ~ 1 (c.f. also Vitello and Salvati, 1976, Goodman, 1986). 
This remains so, however, only until an expansion factor r(p max )/r ~ rj 2 is reached. After 
that, the lab width starts to grow linearly (see Fig. 5a). The comoving frame width, for 
its part, grows initially linearly with the expansion factor until rj, remains approximately 
constant between rj and rj 2 , and then resumes a linear growth with radius (Fig. 5b). The 
physical reason for this behavior is discussed in §3. 
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We have also calculated numerically the lab radius at which an expanding fireball 
becomes Thomson optically thin against its own electrons, as a function of the initial 
radiation energy to rest mass ratio rj. This condition is computed in the comoving frame, 
being the same in the lab frame, as it should, since the optical depth of a constant amount 
of mass is an invariant between reference frames. This (lab) thinning radius r t is found to 
vary as ?? -1 / 2 , in excellent agreement with the analytic estimate (3.13). Another quantity 
of interest is the value of average bulk Lorentz factor at which fireballs of varying rj become 
optically thin, as a function of r\. The numerically computed value of I\ grows linearly with 
rj for low values of rj where optical thinness is achieved before the bulk velocity V saturates 
to ?y, and T t decreases as r? -1 / 2 for larger values of r\ for which thinness is achieved after 
the bulk Lorentz factor has saturated (in agreement with Shemi and Piran, 1990). The 
linear decrease oc n~ 1 / 2 changes to rf for very large values of rj (or initial radiation energies 
sufficiently larger than M a c 2 ) at which the opacity of the fireball at optical thinness is still 
dominated by the pairs (instead of by the polluting baryonic electrons) . This occurs for 
rj ^ 0.6 x 10 10 , above which r t ~ constant, T t ~ 2.4 x 10 3 ~ constant (see §3). 



3.1 Accelerated and Saturated Expansion 

A fireball expanding into a low density external environment may, in its initial stages, 
be considered to be expanding in a vacuum. Thus, as long as the inertia of the accumulated 
external matter can be neglected, shocks will not play an important role in the energetics. 
Pressure gradients may exist in the expanding fireball gas, particularly if the initial density 
and pressure distribution is inhomogeneous, but a free (unimpeded) expansion will tend to 
stretch these out. One may therefore, as a first approximation, neglect internal pressure 
gradients. Thus, even though we are dealing with a fluid, the dynamics of the expansion 
would be expected initially to resemble that of a collisionless gas of relativistic particles. 
The particles initially have an isotropic velocity distribution within a radius r G , with an 
initial thermal Lorentz factor r\ = E /M c 2 . But as they expand, the velocity distribution 
viewed in the lab frame will become increasingly anisotropic. For the ballistic expansion 
in the lab frame, when the particles have reached a radius r their velocity vectors will be 
confined within an angle (r/r G ) _1 of the radial direction. A transformation to a frame of 
reference which is outwardly moving with a Lorentz factor 



makes the particle velocity distribution appear again isotropic. This frame is, therefore, the 
comoving frame of the expanding particles in the initial acceleration phase, which moves 
outward with a velocity increasing according to V oc r. (We have introduced in (3.1) an 
angular factor 9 which approximately accounts for jet geometry; for spherical symmetry 
9 ~ 1, see below). 

The expansion occurs at the expense of the comoving frame thermal energy E, which 
from conservation of energy must therefore drop as E oc r _1 . This agrees with the adiabatic 



3. Analytical Treatment of the Fireball Dynamics 




(3.1) 
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expansion law E oc I/ -1 / 3 oc p 1 ^ 3 , where V, p are comoving volume and comoving baryon 
density, provided that the comoving volume V oc r 3 . The latter follows from the fact 
that the laboratory frame radial extent of the particles in free expansion is expected to 
be (initially, at least) Ar ~ r Q , while the comoving radial width AR is related to the 
corresponding lab width through AR = ArT (we denote comoving radii with capital R 
and laboratory radii with lower case r). Since the dimensions transverse to the motion will 
be the same in the lab and comoving frame, the comoving volume in this (acceleration) 
phase is V oc r 2 AR oc r 2 r Y oc r 3 . 

The accelerating behavior V oc (r/r ), however, can only go as long as the internal 
energy of the expanding gas is relativistic. After the comoving energy density drops below 
the baryon rest mass density pc 2 , the bulk Lorentz factor T must saturate to the maximum 
value it can acquire, which is approximately the initial thermal Lorentz factor rj in the case 
of significant baryon loading, or more generally 

T s ~ min [rj , T m ] , (3.2) 

where the second value is appropriate for rj > Y m given in eq.(3.14), as shown in §3.4. 

The initial acceleration and the saturation behavior can also be obtained from a 
phenomeno logical expression for the average bulk Lorentz factor (e.g., Shemi and Piran, 
1990) 

E + M Q c 2 = 77 + I 
~ E + M c 2 V (E/E ) + 1 

making the analogy with a section of the expanding universe, i.e. assuming spherical 
symmetry, homogeneity and isotropy. Here then E is the fluid thermal energy density 
(initially mainly in photons and leptons, with an admixture of M Q baryons) . As long as it 
is optically thick and rj *S 1, the pressure is radiation-dominated throughout the expansion, 
and using the adiabatic behavior of a relativistic fluid one has 

E/E = T/T = (p/po) 1 ' 3 = (V/Vo)- 1 / 3 , (3.4) 

where E, T, V, p are comoving thermal energy, temperature, volume and baryon density. 
In the linear expansion phase, we can take (from the discussion below equation (3.1)) 
yi/3 ^ r -i^ w h ere r j s radius, so From (3.3) and (3.4) one gets for E a /E <C rj the 
accelerating behavior V ~ E Q /E ~ r/r D , while for E D /E 3> rj one gets the saturated r — > rj, 
c.f. eq.(3.2). The saturation occurs (in spherical symmetry) at r s /r D ~ rj. However, for 
rj > T m given by eq.(3.14), rj must be replaced by T m , see §3.4. 

3.2 Lab and Comoving Geometry 

As the particles move outward with velocity vectors which are increasingly radial, 
they form a radially expanding shell whose radius is initially Ar ~ r D . The radial velocity 
spread (v — c)/c~ T -2 , and a noticeable departure from the approximately constant width 
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r Q starts to become appreciable only after r ^ r b where Ar ~ rAv/c ~ rj,?? 2 ^ r G , or 
r b /r a ^ rf . The laboratory frame width is therefore 

Ar~ m «[r.,r/n~{^ i £ ^ * £ _ (3.5) 
The width in the comoving frame is AR = ArT, which from eqs. (3.1), (3. 2) and (3.5) is 

{02/3 r £ Qr r <^ 

r o r s for r s ^ r ^ r & ; (3.6) 

r/r s for r ^ r b , 

where the two characteristic radii for saturation and for the start of the shell linear expan- 
sion in the lab frame are 

r s /r ~ #- 2 / 3 r s , r b /r ~ r 2 . (3.7) 
The comoving volume in the three different stages of free expansion is 

{4ir9 2 r 3 for r < r s ; 

47r0 2 r s r o r 2 for r s < r < r b ; (3.8) 
47r(9 2 r7 1 r 3 for r > r b ; 

and the comoving radiation energy and temperature vary according to 

r #" 2 / 3 (r /r) for r < r s ; 

( — ) = ( — ) = ^ ^" 2 / 3 (r /r s )(r s /r) 2 / 3 = 9~ 4 / 9 T7 1/3 (rjr) 2 / 3 for r s < r < r b ; 

{ e- 2 / 3 (r /r s )(r s /r b ) 2 / s (r b /r) = Q- 4 / 9 ^ 3 (r /r) for r > r b ; 

(3.9) 

This is therefore the complete r dependence of the adiabatic law (3.4). The factor 9 can 
be taken to include both an angle dependence and a possible dependence on the statistical 
weight factor g in the equilibrium energy density gaT 4 . That is, #" 2 / 3 = (g o /g) 1/3 0'~ 2/3 , 
where (g / 'g) 1 ^ 3 = (11/4) 1 / 3 accounts for the change in the statistical factor for the equilib- 
rium energy density when the pairs annihilate as the temperature drops below kT ^ m e c 2 , 
and 6>'~ 2 / 3 accounts for the possibility of expansion along a restricted range of solid an- 
gles. If we have radial expansion along two jets of solid angle O , and all of the energy is 
channeled into this solid angle, then 6' ~ 9 2~ 1 / 2 (or ~ o /2 if there is one jet only). If 
the gas uses some its internal energy to perform work against a medium which prevents 
escape in the directions outside of the jet angles 6* , then the effective 9' < O 2~ 1 / 2 . In 
the discussion below (and in the numerical calculations of §2) we have ignored the small 
shift of the curves associated with the change of the statistical factor, i.e. we simplify to 
{9o/g) ~ 1, 9 ~ 9' . (Note that this angular factor is only approximately correct, since it 
does not account for any possible lateral expansion or transverse radiation loss of the jet; 
it is strictly valid for one-dimensional radial variations within the angle 9'). For a jet, the 
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asymptotic accelerated behavior has the same radial dependence, but displaced to higher 
radii by 6> _2 / 3 , as is the saturation radius, 

T~e 2/3 (r/r ) , r s /r ~6- 2/3 V , T s ~ V , . (3.10) 

(which is valid only for 77 < T m ; otherwise r\ must be replaced by T m given by eq.(3.14) 
see §3.2). The effect of the approximate angular factor, e.g. in the behavior of the growth 
of the bulk Lorentz factor (eq.(3.1)), is that, for the same initial energies the jet achieves 
the same velocities as a corresponding spherical flow, but at larger radii. This is because 
at the same radii the jet, being confined to narrower angles, has a larger internal entropy 
than the spherical flow. The saturation Lorentz factor is however the same in both cases. 



3.3 Optical Depth and Saturation Lorentz Factor 



Depending on the value of n = E /M c 2 , the scattering depth of the fireball at 
optical thiness is dominated by the "polluting" baryonic electrons (if n < n p , where the 
latter is given by eq.(3.19)), or, if the fireball has very low baryon pollution, the scattering 
opacity previous to reaching optical thinness is dominated by the e pairs (Goodman, 
1986, Paczynski , 1986, Shemi and Piran, 1990), for n > rj p . 



a) Electron-dominated scattering: In the case r\ < rj p , when baryonic electrons dominate 
at thiness, the optical depth is 

^ = M /m p aAr = , E k ,,r , 1 = g _ 4/3 r^ ,r a . 2 ^ n . 

47rr 2 9 2 r 2 Ar 4ivr 2 c 2 r 9 2 n r\ r 

Here a is scattering cross section, k = a/m p ~ 0.4 cm 2 g -1 , and T m is defined in eq.(3.14). 
The initial optical depth is r G = r^ -4 / 3 ?? -1 = T, k = 77 _1 E r)0 K, where 

E = M /4nr 2 9 2 , E r , = E c~ 2 /4nr 2 a 9 2 , (3.12) 

are the initial baryon mass surface density and the initial radiation equivalent mass surface 
density. In the course of the expansion, the gas becomes optically thin at a radius r t defined 
by 

I* = 9- 3 / 2 rU 2 V~ 1/2 = (S «) 1/2 = r l J 2 = 1.9 x lO^r^-V 172 , (3.13a) 

or 

_ , E k a/2 _ 
Tt ~ ^A^c 2 ^ 

where the initial energy E Q and radius r Q have been arbitrarily normalized to 10 51 ergs and 
10 6 cm, and the second version (3.13b) is independent of r a . The critical Lorentz factor 
T m and the critical rj m are defined as 

r m = Vm = 9^(r r ] ) 1 / 3 = 9^{Y, nr 1 )^ = 9^{^ k)^ 

= (_^l_) 1 /3^-2/9 = 3 3 x 1Q 5 E V3 r -2/3^-2/9 ^.lA) 

QTtr a c 
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This is the maximum possible bulk Lorentz factor achievable for a given initial radiation 
energy E Q deposited within a given initial radius r Q (§3.4). It is also, for a given E and r Q , 
the "critical" value of r\ = rj m (reached at a critical loading mass M Q = M m = (E /rj m c 2 )) 
for which the thinning radius (3.13) is equal to the saturation radius (r s /r Q ) ~ 9~ 2 / 3 rj. For 
a given E Q , as one varies M Q or ?y, the critical value r\ = rjm = T m is reached at a radius 
r m given by 

t jil = a = a „ 9-^r m . (3.i5) 

r a r a r a 



For rj < r m , optical thinness is reached after saturation, whereas for r\ > Y m optical thin- 
ness is reached before saturation, so the bulk Lorentz factor at thinning T t = 2 ^ s (r t /r o ) 
is equal to 

r t = [ v , r^ /2 ?r 1/2 ] for [ v < r m , r m < v < r p ] , (3.16) 

where T p is given by eq.(3.20). 

b) Pair-dominated Scattering: For very large rj, such that rj > 1] P > T m , the e^ 1 scattering 
dominates, instead of baryonic electron scattering, and optical thiness occurs when the 
pairs fall out of equilibrium, which occurs at a lab radius r p given by 

r JL = fl- 2 ^ „ 2.4 x 10 3 El{\- 3/4 e- 2 / 3 , (3.17) 

where T p ~ 15 keV is the temperature where the thermal pair density becomes negli- 
gible, and T a ~ 36.4i£g{ 4 r 6 3 ^ 4 MeV is the initial temperature. This occurs during the 
acceleration stage, so at this point, where pairs become optically thin, V has reached the 
value 

T t = T p = E /E p = T /T p = e 2 '\r p /r ) ~ 2.4 x ltfE\{ A rf* , (3.18) 

valid for rj > rj p . The latter is the value of r\ above which the pair-dominated regime 
occurs, 

V P = rl/T 2 p ~ 0.63 x lO 10 ^" 1/2 »- 2/3 , (3-19) 

in terms of which one can write the Lorentz factor and the radius at which pairs become 
optically thin as 

F 3/2 r 3/2 

V p = % , r f = 0-^Y p = 0- 2 ^. (3.20) 
Vp r ° V P 



3.4 Photon Drag and Final Baryon Lorentz Factor 

Photons are obviously coupled to the baryons when r > 1 , which ensures a radiation- 
like equation of state with adiabatic index 4/3. However, even after r < 1, baryons may 
remain coupled to the photons if the density of the latter is so large that the Compton drag 
time is shorter than the expansion time in the comoving frame. The comoving Compton 
drag time is the time during which an electron sweeps up an amount of photons whose 
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mass equivalent equals one proton plus one electron's mass, to ~ m p c 2 / arcu 1: where -u 7 
is the comoving photon energy density. This is longer than the Compton cooling time by 
a factor m p /m e . The ratio of the comoving baryon rest mass energy density u m and the 
comoving photon energy density w 7 is 

£ = S = ^? = , f' (3 ' 21) 

which is equal to 6~ 2 / 3 rjr /r if r < r s . Dividing the Compton drag time to by the 
comoving expansion time AR/c and multiplying and dividing by the comoving baryon 
density n (comoving baryon rest energy density u m ), the ratio of the Compton drag time 
to the expansion time is 

^^h^¥^vf^^vf' <3 ' 22) 

where r is the optical depth (3.11), and the last equality is valid for r ^ r s , with different 
dependences for r ^ r s . The final decoupling of photons and baryons occurs at a radius 
rf where £ = 1, or at the thinning radius rt, whichever is largest. 

The final bulk Lorentz factor of the baryons is the value that it has at decoupling 
from the photons, either when the fireball becomes optically thin or when the Compton 
drag has become longer than the expansion time, ( > 1, whichever occurs last. For rj < T m , 
the fireball saturates to a value r = 77, and thinning occurs at a radius larger by a factor 
(ft/f s ) = (Tm/v) 3 ^ 2 > 1 where Q = (T m /v) > 1) so the final baryon Lorentz factor is 
Tf = rj in this case. On the other hand, for rj > T m , at optical thinness St = (??/r m ) 3 / 2 > 1 
and Qt = (^m/v) 3 ^ 2 < 1? so the Compton drag keeps photons and baryons coupled beyond 
r t . This means that, even though most photons do not collide, most electrons do keep 
colliding with a small fraction of the photons, and keep being accelerated (with their 
baryons). They finally decouple when Q = 1 at (r//r D ) = 9~ 2 ^ 3 T m , where Tf ~ T m , 
Ef = (rj/Tm) > 1 and r/ = (T m /rj) < 1. The terminal baryon bulk Lorentz factor is 
therefore 

T f = min [ V , r m ] , for [77 < T m , V > T m ] . (3.23) 

This is equal to the Lorentz factor at thinning T t in the case rj < T m , when T t = Tf = rj, 
but Tt 7^ Tf for rj > T m . 

3.5 Baryon Loading and Final Radiation to Kinetic Energy Ratio 

One can distinguish four baryon-loading regimes, characterized by the value of n: 
H) High-load fireballs , for 1< 77 < T m , or 0.6 x 1O" 3 M £ M a £ 1.8 x 1O" 9 M . In this 
regime T grows linearly with r/r a until reaching r s /r a = 9~ 2 ^ s rj, where it saturates to n, 
and becomes optically thin at r t /r a = r m (r m /?7) 1 / 2 > T m > r s /r a . The observer-frame 
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(Doppler blueshifted) energy of the radiation escaping under the assumption of isotropy 
(even if in reality it is beamed) and the kinetic energy of the baryons is 

E r)0b = E 9~ 4/9 r]~ 1/3 ( — ) 2/3 = E (r^—) < E , E k = M c 2 r, ~ E Q . (3.24) 

That is, the thinning radiation burst is a small fraction of E Q (unless r\ — > T m ), most of 
the energy having gone into the kinetic energy of the baryons. 

M) Critical-load Fireballs, for rj = T m or M Q ~ 1.8 x 10 -9 Mq. In this case T s = T m at 
r s = rt = ff, and both the observer- frame radiation energy and the baryon kinetic energy 
assume their maximal value, E rjQ b ~ E Q /2, Ek ~ E Q /2. 

L) Low-load or underloaded fireballs, for T m < rj < rj p , or 1.8 x 1O _9 M & M Q ^ 0.95 x 
10 -13 Mq. For these (rather low) values of the baryon pollutant mass, V grows linearly 
with r/r Q until becoming optically thin at rt/r Q < r s /r Q ~ r m /r , where it has the value 
r t ~ 6 2 / 3 r t /r = r m (r m /?7) 1 / 2 < rj. Most of the radiation energy escapes from here 
without further scattering, while the baryons continue coupled by Compton drag to a small 
fraction of the photons until the radius r/ where Ff = T m . The bulk of the observer- frame 
radiation energy observed (assumed over 4%) and the final kinetic energy of the protons 
are 

E r , ob ~ E 6- 2 / 3 (^)T t ~E , E k = M c 2 T f ~ £ c ( — ) < E D . (3.25) 

P) Pair fireballs, for r\ > r\ v or M Q ^ 0.95 x 1O _13 M . In this extremely underloaded 
regime, the pairs dominate the Thomson opacity and the value of V grows linearly with 
r/r Q but it always becomes optically thin at the same value of the radius where r ~ T p , 
r t /r D = r p /r a ~ 6~ 2 ^ 3 T pi given by eqs.(3.18)-(3.20), independent of rj as long as the latter 
is greater than rj p . The observer-frame radiation energy at thiness under the assumption 
of isotropy and the final kinetic energy of the baryons is again given by eq.(3.25). 
The significance of the ratio of observed radiation energy to kinetic energy, in our "stan- 
dard" model (Meszaros and Rees, 1993), is that this represents the ratio of the energy in 
a short precursor burst to that in the main burst, the latter coming from the recovery of 
the kinetic energy upon interaction with the external medium (§4). 

4. Magnetic Fields, Efficiency and Photon Energies 

4-1 Magnetic Fireballs 

In almost any model of the initial energy release, the initial mass motions are expected 
to be extremely violent (v ~ c), and this could magnify any pre-existing magnetic fields, 
via compression, shearing, turbulent dynamo mechanisms, Parker-type instabilities, etc. 
(e.g. Usov, 1992, Narayan, Paczyhski and Piran, 1992, Thompson and Duncan, 1993). If 
the initial total "disposable" energy (i.e., that portion of the available gravitational energy 
that is not lost in the form of neutrinos) is E a = 10 51 i?5i ergs, this might be distributed 
between radiation and magnetic components as 



Ebo — £,E , E ro — (1 — O-^o , 



(4.1) 
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where £ < 1. One may define separate magnetic and radiation rj parameters 

VB = E Bo /M c 2 = , Vr = E ro /M Q c 2 = (1 - OV (4-2) 

in terms of the usual total t] = E /M c 2 , with 1] = 1] B + rj r = £77 + (1 — Even 
if the magnetic fields have a large scale structure, the expansion is essentially isotropic in 
the comoving frame, and therefore the magnetic energy density will evolve oc y -4 / 3 , so 
the total magnetic energy and radiation energy in the comoving frame vary in a similar 
manner, 

Eb/Ebo = B 2 /B 2 Q = E r /E ro = (V/K)" 1/3 , where the volume factors are given 
by eqs. (3.9). The bulk Lorentz factor, similarly to (3.3), is again 

r ^ E + M 7/ + 1 77+1 , s 

E r + E B +M a rj r (E r /E ro ) + ri B (E B /E Bo ) + 1 e- 2 /%r /r)rj + 1 ' 1 ' } 

so that, as in the case of pure radiation, V ~ 6> 2 / 3 (r/r ) for r & r s while V ~ Tj for v ^ T S j 
with r s /r Q = 6~ 2 ^ 3 T S and T s = min[?y, T m ], where T m is given by eqs(3.14). The behavior 
is therefore the same as in §3, even if all the disposable energy E goes into magnetic fields 
(£ — ► 1), i.e. a purely magnetic fireball. The maximal initial magnetic field will be 

B = (8ntE /V o y/ 2 = l^e /2 El{ 2 rf /2 G , (4.4) 

and the comoving field strength evolves according to 

C e-^{r Q /r) 2 for r < r s ; 

(^-) = { 9- 8 / 9 T7 2/3 (r /r)^ for r s < r < r b ; (4.5) 
{6-^T 2/3 (r /r) 2 for r > r & ; 

where r s , r& are given in (3.7). Even after most of the pairs annihilate, the exponentially 
small fraction of frozen-in surviving pairs is enough to provide the currents needed to 
support large-scale fields, so that the behavior (4.5) is uninterrupted beyond the radius at 
which annihilation, optical thinness, photon drag decoupling, etc. occur. 

We may also consider briefly the most extreme scenario of magnetic field dominance, 
that where the initial field B Q is in equipartition not with the disposable energy E ~ 
10 51 -E5i ergs (which, as in supernovae, is of order 10~ 3 of the total liberated binding 
energy from the gravitational collapse of about a solar mass) but rather is in equipartition 
with the binding energy itself, E h ~ GMq/R n ~ 10 54 E 54 ergs. This might occur in 
a rapidly-spinning object, where rotational and shearing motions would involve a large 
fraction of the total binding energy. In this case 

B = B m - 3 x 10 18 (£/30) G, (4.6) 

where, in terms of the first of equations (4.1), £ ~ 30. Now, however, r\ r is no longer 
defined as (1 — £)?y but rather just as rj r = E ro /M c 2 . The dynamic considerations are 
similar to those just described, but the dynamics is entirely described by using t]b instead 
of i] everywhere in §3, and increasing T m in eq.(3.14) by a factor 10 (since E 5 i ~ 10 3 ). 
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4-2 Standard Shock Deceleration Model 

The magnetic fields are carried on in the comoving frame of whatever polluting 
baryons M D were present originally in the fireball, and continue expanding with the latter 
and the corresponding electrons plus surviving pairs until the fireball matter is decelerated 
by the external medium (Rees and Meszaros, 1992). For a fireball which has saturated its 
bulk Lorentz factor to a value T s = min[r/, T m ] the deceleration occurs when the mass of 
external matter swept up by the fireball equals T~ 1 M ; the swept-up external matter is 
shock heated to a comoving average thermal Lorentz factor 7 S ~ F s ~ rj » 1, while a 
reverse shock starts to propagate into the adiabatically cooled fireball material, eventually 
reheating it to a marginally relativistic average thermal Lorentz factor 7 r ~ T r & 2. If the 
forward and reverse shocks are able to create a field which approaches equipartition with 
the post-shock thermal energies of the particles (similarly to what appears to occur in su- 
pernova remnant shocks), both the forward and reverse shock regions can radiate away the 
entire thermalized bulk kinetic energy of the shock, i.e. the entire total initial disposable 
energy E Q of the event. However, because of the dependence of the radiative efficiency on 
the magnetic field strength, it is worthwhile to consider departures from shock equiparti- 
tion, as well as the effect of frozen-in primordial fields discussed in the previous subsection, 
whose initial strength B Q may have been amplified by shear or turbulent dynamo effects 
during the cataclysmic event of the initial energy release, parametrized through £ defined 
in the first of equations (4.1). 

In order to estimate the radiative efficiency at the deceleration shock, we use the 
parameters of our "standard" shock model, e.g. Meszaros and Rees, (1993). For typical 
parameters, the shock deceleration occurs after the fireball has saturated (and after it 
has become optically thin, producing a brief and weak precursor burst). The deceleration 
radius at which the fireball gas starts to feel the inertia of an external medium of 
uniform density n ext = n x cm" 3 is r d /r Q = (r/ ex t/V) 1/3 , where rj ext = (E /V n ext m p c 2 ), 
or ra = (3E /Ann ext m p c 2 ri 2 ) 1 ^ s 1 which is 

ml6 1/3 ,-,1/3 -2/3 / * n\ 

r d ~ 10 iO V E 5 [ 7/3 1 cm . (4.7) 

(For maximal magnetic dominance, E 5 i ~ 10 3 , rj = tjb , £ ~ 30, but for the purposes 
of this example we consider the more conservative case of E 5 i ~ 1 , £ ^ 1 , n\ ~ 1). 
The comoving width of the forward and reverse shocked shells is r\ times larger than their 
laboratory width ~ r^jn 2 , or AR ~ ra/n ~ 10 13 r/ 3 d ^ 3 cm. The comoving expansion time 
is consequently 

t ex = At c = AR/c ~ lO 3 ^ 573 s , (4.8) 

while the laboratory expansion time (the observed GRB burst time, if the shock is radia- 
tively efficient) is At = ra/n 2 ~ 1 r/ 3 8 ^ 3 s. The total number of baryons (protons) involved 
in the fireball is N PjQ = E /nm p c 2 ~ 10 54 ?y _1 , and the comoving density in the preshocked 
fireball gas (for > r& in eqs.(3.8) with (3.4)) is 



n = n d = N Pi0 n/Aixe 2 rl ~ lO 6 ^ cm" 3 



(4.9) 
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If the reverse shock, whose bulk Lorentz factor achieves at most a value F r ~ 2, produces a 
shock equipartition field, this is Bd e = (8irnm p c IV4) 1 / 2 - 10 2 ?73 G, while if the magnetic 
field is mainly the original, possibly amplified and adiabatically expanded fireball magnetic 
field (4.4) or (4.6), this is 

B d = AB r, 2 /\r d /r ) 2 = 4 x lO" 1 ^ 2 G , (4.10) 

where in both cases we have included a factor 4 to take into account the shock com- 
pression factor. Aside from the different rj dependence, we may consider both the shock- 
equipartition and frozen-in fields to be given by (4.10), where the shock-equipartition field 
case has £ ~ 250?y^" 1 while the frozen-in field has £ & 1 (or £ ^ 30 in the extreme magnetic 
dominance case). 

4-3 Radiative Efficiency and Photon Energies 

The efficiency of the blast wave moving ahead of the contact discontinuity can only 
be affected by the field B d e developed in the shock, as discussed in our previous papers. 
Here we focus on the radiative efficiency of the reverse-shocked fireball material, which can 
be affected by the original, possibly amplified, fireball field Bd of eq.(4.10). The "fireball" 
at the stage just before the reverse shock moves into it contains only a small fraction of 
the cooled original pairs and the original polluting baryons and electrons. With a field 
of order (4.10), relativistic electrons of Lorentz factor 7 ~ 10 6 produced by the reverse 
shock are sufficient to ensure a synchrotron cooling time comparable to the (comoving) 
expansion time (4.6), i.e. near unit radiative efficiency, and photons of characteristic 
energy ~ I773 MeV in the laboratory frame. If shock acceleration is considered, limited by 
synchrotron losses, the maximum possible 7 ^ 2 x 10 7 B~ 1 / 2 and the maximum energy of 
the synchrotron photons is ~ 1773 GeV in the laboratory frame, independent of the field 
strength. 

Synchrotron losses, however, will almost certainly be surpassed in importance by 
the inverse Compton (IC) losses of the electrons, in the radiation field of the synchrotron 
photons. Without diffusive Fermi-type acceleration , the electrons passing through the 
shock may energetically attain a maximum average Lorentz factor 

7 ~ (m p /m e )T r - 4 x 10 3 C . (4.11) 

With this average Lorentz factor, the ratio of synchrotron cooling to expansion time 
tsy/tex ~ 5 x 10 8 7 _1 S~ 2 ~ 2.5 x 10 3 C _1 £~ 2 ?73~ 7 / 3 ^ 1> i- e - the synchrotron efficiency 

1 /3 

is low. The energy density in synchrotron photons is u sy ~ nP sy AR/c ~2x lO 2 ^ 2 ^ ub, 
where ub = B 2 /87T is the magnetic energy density. The ratio of inverse Compton cooling 
time to expansion time is therefore 
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giving a first-order inverse Compton efficiency of 10% even for this modest average Lorentz 
factor 7 ~ 4 x 10 s ( and £ ~ 1. The synchrotron photons (IR in the comoving frame, optical 
in the lab frame for this 7 and 773 ~ 1) are boosted up by the single inverse Compton 
scattering to MeV and GeV energies in the comoving and lab frames. 

The radiative efficiency will be larger in the presence of a diffusive (or similar) shock 
acceleration mechanism that produces an electron power-law distribution. Considering the 
electron shock acceleration to be limited by first-order inverse Compton losses (using the 
Thomson limit in the electron rest frame and balancing tic with the acceleration time 
t a ~ 1 I^b where ujb = eB/m e c) one gets a maximum 7 significantly higher than (4.11) 
by a factor £ m ~ 4.5 x lO 1 ^ -1 / 4 ^ 7 / 12 . The synchrotron photon energy in the electron 

rest frame becomes of order 1 MeV at a slightly lower energy, ( m ~ 2 x lO 1 ^ -1 / 3 ^ 2 ^ 3 . 
Since the transition to the Klein-Nishina cross section at this energy will effectively limit 
the maximum Lorentz factor to this value, we will use this second value of £ m ~ 20, i.e. 
7 m ~ 10 5 . With this, the synchrotron and first-order IC photons have a characteristic 
maximum energy of 

E sy ~ 10e 1/3 r?3 /3 k eV , E IC ~ 10 14 r 1/3 % /3 eV , (4.13) 
in the lab frame, and a radiation efficiency 

e sy ~ 10- 2 £ 5 / 3 77 3 5/3 , e IC ~ min[l , lO 3 ^ 3 ] . (4.14) 

Another effect which could increase the efficiency is the higher order IC scattering 
(e.g. Rees, 1967). The Thomson optical depth of the fireball shell material is 

t t ~ na T fAR ~ 0.6 x 10" 5 /% /3 , (4.15) 

where / < 1 is the fraction of the fireball material that has been heated by the shock. 
When tt is larger than 7 -2 (as long as the scattering is not Klein-Nishina dominated) the 
conditions are fulfilled for both first order and higher order inverse Compton scattering to 
dominate over synchrotron losses. This is because the synchrotron, first-order IC, second- 



order IC, etc. photon energy densities are 

u sy ~ nP S yfAR/c ~ na T ^ 2 u B fAR ~ t t ^ 2 u b , (4.16a) 

uic ~ nPicfARjc ~ na T ^ 2 u sy fAR ~ t^ 4 u b , (4.166) 

uic ~ nPjcfAR/c ~ naT^ 2 ujc f AR ~ t^^ub , ••• (4.16c) 



and for tt ^ 7 -2 the ratio of the successive energy densities in eq(4.16) are larger by 
increasing factors of r^7 2 > 1- Even for 7 ~ 10 3 the IC losses exceed synchrotron losses by 
a factor 10, and the second order IC losses would exceed first order IC losses by another 
factor of 10, were it not for the fact that for second order IC and the maximal 7 the photons 
in the electron rest-frame are now subject to the Klein-Nishina decrease in the cross section. 
(For 7 ~ 10 5 the respective factors are 10 5 , 10 10 ,, etc., and the importance of the IC losses 
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make themselves felt already when the shock has heated a small fraction / *5 10~ 5 of 
the fireball material). Notice that for both our simplest model (using the average Lorentz 
factor 7 ~ 10 3 ) and for the non-thermal (diffusive acceleration) model with r ) m ~ 10 5 , the 
first order IC is in the Thomson limit and exceeds synchrotron losses, but the Klein-Nishina 
losses diminish the importance (at least as far as the efficiency is concerned) of the second 
and higher order IC scattering, although they will affect the spectrum. In the power-law 
spectrum model, the first-order IC efficiency is unity for > 10 -3 , i.e. for initial fields 

B Q i<S 10 14 ?73— 2/3 G in the fireball, or for fields developed by turbulent instabilities in the 

_ 2/3 

reverse shock which are a factor 10~ rj 3 below equipartition with the thermal (reverse) 
post-shock particles, a fairly undemanding assumption 

5. Discussion 

We have presented both analytical and numerical calculations of the evolution of 
a freely expanding ultrarelativistic gas produced by an impulsive energy release. The 
expansion and cooling from an initially optically thick and extremely hot fireball is followed 
through the stage where the rest- frame energy density becomes comparable to the rest mass 
density, where the expansion bulk Lorentz factor, which until then grew proportionally to 
the expansion factor, saturates to a value which is equal to the smallest of the ratio of 
the initial thermal energy to rest mass energy in the fireball or the bulk Lorentz factor at 
which Compton drag becomes negligible. The laboratory and comoving frame geometry 
as well as the radiation and gasdynamic variables were investigated for a large range of 
values of r\ = E /M c 2 , and we derived analytic scaling laws showing the dependence of the 
variables as a function of the various initial parameters, including an approximate angular 
scaling for the case of jet-like expansion within a limited range of solid angles. 

We have also devoted particular attention to the likely magnetic field content of fire- 
balls. An initial magnetic field, possibly amplified by the violent mass motions in the initial 
impulsive event, may contribute a dynamically-significant fraction of the total fireball en- 
ergy; there is also the possibility that a field develops in the shock heated fireball after 
the ram pressure of an external medium decelerates the expansion. When the fireball runs 
into an external medium and is decelerated, the efficiency with which the re-randomized 
energy is radiated depends on the magnetic field strength, and on how the electrons are 
accelerated by the resultant shocks. (The radiation processes during the deceleration phase 
which give rise to the typical burst profiles are, fortunately, rather less sensitive to precisely 
how the fireball was originally formed - the formation mechanism may have resembled the 
impulsive model discussed in sections 2 and 3, or could alternatively have been spread over 
as long as a second, as in the scenarios of Usov (1992) or Narayan et al (1992) ). We find, 
for a range of assumptions, that inverse synchro-Compton cooling should be extremely 
efficient in radiating away the kinetic energy of the cooled fireball in a time scale short 
compared to the expansion time, which gives the right order of magnitude duration and 
total energy for a gamma ray burst source, and produces a nonthermal spectrum with the 
bulk of the energy in the the MeV to GeV range. 
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Figure Captions 

Fig. 1 : Mass distribution in the laboratory frame for a fireball of r\ = 10 4 at the times 
when the lab expansion factor r(p max )/r ( proportional to the lab time) is equal to 
0.2, 0.4, 0.6, 0.8, 1.0 times rj, from top to bottom. 

Fig. 2 : Distribution of the bulk Lorentz factor F over the various mass shells (i.e. 
against the Lagrangian mass coordinate), for r\ = 10 , at the times when the 
lab expansion factor ( which is proportional to the lab time) is r(p max ) jr = 
r] 1 / 2 , 7], r? 3 / 2 , rf, if/ 2 \ if, from bottom to top. 

Fig. 3 : The average bulk Lorentz factor T as a function of the expansion factor r(p max ) / r D , 
for r] = 10 4 . 

Fig. 4 '■ a ) The lab frame mass density against the Lagrangian mass coordinate, 
at the times (from top to bottom) when the expansion factor is r(/9 max )/r = 
77 1 / 2 , ?7, if/ 2 , rj 2 , if/ 2 , if. 

b) The comoving frame mass density against the Lagrangian mass coordinate, for 
the same instants as a) 

Fig. 5 : a) The lab frame width of the shell where most (80%) of the mass is located, as 
a function of the expansion factor r(p max )/r , for r\ = 10 4 . 

b) The comoving frame width of the shell containing most (80%) of the matter, as a 
function of the expansion factor r(p max )/r , for rj = 10 4 . 



